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Can Inflating Braneworlds be Stabilized? 
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We investigate scalar perturbations from inflation in braneworld cosmologies with extra dimen- 
sions. For this we calculate scalar metric fluctuations around five dimensional warped geometry 
with four dimensional de Sitter slices. The background metric is determined self-consistently by the 
(arbitrary) bulk scalar field potential, supplemented by the boundary conditions at both orbifold 
branes. Assuming that the inflating branes are stabilized (by the brane scalar field potentials), 
we estimate the lowest eigenvalue of the scalar fluctuations - the radion mass. In the limit of 
flat branes, we reproduce well known estimates of the positive radion mass for stabilized branes. 
Surprisingly, however, we found that for de Sitter (inflating) branes the square of the radion mass 
is typically negative, which leads to a strong tachyonic instability. Thus, parameters of stabilized 
inflating braneworlds must be constrained to avoid this tachyonic instability. Instability of "stabi- 
lized" de Sitter branes is confirmed by the BraneCode numerical calculations in the accompanying 
paper pj. If the model's parameters are such that the radion mass is smaller than the Hubble 
parameter, we encounter a new mechanism of generation of primordial scalar fluctuations, which 
have a scale free spectrum and acceptable amplitude. 



I. INTRODUCTION 

One of the most interesting recent developments in 
high energy physics has been the picture of braneworlds. 
Higher dimensional formulations of braneworld models in 
superstring/M theory, supergravity and phenomenologi- 
cal models of the mass hierarchy have the most obvious 
relevance to cosmology. In application to the very early 
universe this leads to braneworld cosmology, where our 
3+1 dimensional universe is a 3d curved brane embedded 
in a higher-dimensional bulk @. Early universe inflation 
in this picture corresponds to 3+1 (quasi) de Sitter brane 
geometry, so that the background geometry is simply de- 
scribed by the five dimensional warped metric with four 
dimensional de Sitter slices 



ds 2 = a 2 (w) [dw 2 ~ dt 2 + e 2Ht dx 2 ] 



(1) 



For simplicity we use spatially flat slicing of the de Sitter 
metric ds 2 . The conformal warp factor a(w) is deter- 
mined self-consistently by the five-dimensional Einstein 
equations, supplemented by the boundary conditions at 
two orbifold branes. We assume the presence of a sin- 
gle bulk scalar field ip with the potential V(ip) and self- 
interaction potentials U±(ip) at the branes. The poten- 
tials can be pretty much arbitrary as long as the phe- 
nomenology of the braneworld is acceptable. The class of 
metrics Q with bulk scalars and two orbifold branes cov- 
ers many interesting braneworld scenarios including the 
Hofava-Witten theory 0,13 i the Randall-Sundrum model 
[UliJ with phenomenological stabilization of branes 0,0, 
supergravity with domain walls, and others 0,^3- 



We will consider models where by the choice of the 
bulk/brane potentials the inter-brane separation (the so- 
called radion) can be fixed, i.e. models in which branes 
could in principle be stabilized. The theory of scalar 
fluctuations around flat stabilized branes, involving bulk 
scalar field fluctuations Sip, scalar 5d metric fluctuations 
and brane displacements, is well understood [ll|- Simi- 
lar to Kaluza-Klein (KK) theories, the extra-dimensional 
dependence can be separated out, and the problem is re- 
duced to finding the eigenvalues of a second-order differ- 
ential equation for the extra-dimensional ( to-dependent ) 
part of the fluctuation eigenfunctions subject to the 
boundary conditions at the branes. The lowest eigen- 
value corresponds to the radion mass, which is positive 
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fluctuations around flat stabilized branes are also stable. 

Brane inflation, like all inflationary models, generates 
long wavelength cosmological perturbations from the vac- 
uum fluctuations of all light (i.e. with mass less than the 
Hubble parameter H) degrees of freedom. The theory of 
metric fluctuations around the background geometry (JTJ 
with inflating (de Sitter) branes is more complicated than 
that for the flat branes. For tensor fluctuations (gravita- 
tional waves), the lowest eigenvalue of the extra dimen- 
sional part of the tensor eigenfunction is zero, m = 0, 
which corresponds to the usual 4d graviton. As it was 
shown in |2ll |22| i massive KK gravitons have a gap in 
the spectrum; the universal lower bound on the mass is 

m > y^H. This means that massive KK tensor modes 

are not generated from brane inflation. Massless scalar 
and vector projections of the bulk gravitons are absent, 
so only the massless 4d tensor mode is generated. 

Scalar cosmological fluctuations from inflation in the 
braneworld set ting (fill have been considered in many im- 
portant works |ll |3 E E [H E 13 H3 Th e theory 
of scalar perturbations in braneworld inflation with bulk 
scalars is even more complicated than for tensor pertur- 
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bations. This is because one has to consider 5d scalar 
metric fluctuations and brane displacements induced not 
only by the bulk scalar field fluctuations Sip, but also by 
the fluctuations 8\ of the inflaton scalar field x living at 
the brane. In fact, most papers on scalar perturbations 
from brane inflation concentrated mainly on the inflaton 
fluctuations Sx, while the bulk scalar fluctuations were 
not included. This was partly because in the earlier pa- 
pers on brane inflation people considered a single brane 
embedded in an AdS background without a bulk scalar 
field, and partly because for brancworlds with two sta- 
bilized brancs there was an expectation that the fluctua- 
tions of the bulk scalar would be massive and thus would 
not be excited during inflation. 

In this letter we focus on the bulk scalar field fluctu- 
ations, assuming for the sake of simplicity that the in- 
flaton fluctuations Sx are subdominant. We consider a 
relatively simple problem of scalar fluctuations around 
curved (de Sitter) branes, involving only bulk scalar field 
fluctuations Sip. We find the extra-dimensional eigenval- 
ues of the scalar fluctuations subject to boundary con- 
ditions at the branes, focusing especially on the radion 
mass m 2 for the inflating branes. In particular, we inves- 
tigate the presence or absence of a gap in the KK spec- 
trum of scalar fluctuations in view of the tensor mode 
result. Our results are a generalization of the known re- 
sults for flat stabilized branes which we reproduce 
in the limit where the branes are flattening H — > 0. 

II. BULK EQUATIONS 

The five-dimensional braneworld models with a scalar 
field in the bulk are described by the action 

S = Ml J ^d 5 x {R- {Vip) 2 -2V{ip)} 

-2MIJ2 J V^d\ {[AC] + £%)} , (2) 

where the first term corresponds to the bulk and the sum 
contains contributions from each brane. The jump of the 
extrinsic curvature [JC] provides the junction conditions 
across the branes (see equation (|llf) below) . Variation of 
this action gives the bulk Einstein Gab = Tab{w) an d 
scalar field Dip — V v equations. For the (stationary) 
warped geometry Q they are 

ip" + 3— cp' - a 2 V = 0, (3a) 
a 

T-'S-^-T- (3b » 

•(£-*■)-£-•*. <*> 

where the prime denotes the derivative with respect to 
the extra dimension coordinate w. The first two equa- 
tions are dynamical, and the last is a constraint. The 



solutions of equations were investigated in detail in 

M- 

Now we consider scalar fluctuations around the back- 
ground Q). The perturbed metric can be written in the 
longitudinal gauge as 

ds 2 = a(w) 2 [(1 + 2<P)dw 2 + (1 + 2^)dsl] . (4) 

The linearized bulk Einstein equations and scalar field 
equation relate two gravitational potentials &(x A ), 
^(x A ) and bulk scalar field fluctuations Sip(x A ). The 
off-diagonal Einstein equations require that 

* = "2' (5) 

similar to four-dimensional cosmology, although the co- 
efficient is different. 

The symmetry of the background guarantees separa- 
tion of variables, so that perturbations can be decom- 
posed with respect to four-dimensional scalar harmonics, 
e.g. 

$(/)=^f m HQ m ((,x), (6) 

m 

where the eigenvalues m (constant of separation) appear 
as the four-dimensional masses A DQ m = m 2 Q m , where 
4 D is the D'Alembert operator on the 4d de Sitter slice. 
The four-dimensional massive scalar harmonics Q m can 
be further decomposed as Q m (t,x) = J f^\t) e l ^ x d z k. 
The temporal mode functions (i) obey the equation 

/ + 3Hf + (e- 2Ht k 2 + m 2 )f = 0, (7) 

where dot denotes time derivative, and we dropped the 
labels k and m for brevity. 

Out of the remaining linearized Einstein equation we 
get the following equations for the extra-dimensional 
mode functions <& m (w) and Sip m (w) 

(a 2 ®)' = H V^, (8a) 

where we again omitted the label m for transparency. 

These are very similar to the scalar perturbation equa- 
tions in four-dimensional cosmology with a scalar field 
|23| . except for some numerical coefficients and powers 
of a(w) (because the spacetime dimensionality is higher), 
and up to time to extra spatial dimension exchange. In- 
deed, we can introduce the higher-dimensional analog of 
the Mukhanov's variable. However, in the presence of 
the curvature term H 2 , the eigenvalue m 2 enters the sec- 
ond order equation for it in a complicated way, similar 
to that in the 4d problem with non-zero spatial curva- 
ture, see e.g. [24j. We can introduce another convenient 
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variable u r . 



3/2 



^ > 



Then the two first order 



2 ip' 

differential equations (JSJ) can be combined into a single 
Schrodinger-type equation 



AH 



V e s(w) ) U m = 



(9) 



|<// 2 where we 



with the effective potential V c s = - 

defined z = (^aip' 2 ) 1 ■ 

There are two main differences relative to the four di- 
mensional cosmology. First, in the latter case, FRW 
geometry with flat 3d spatial slices is usually consid- 
ered, while the five dimensional brane inflation metric 
has curved 4d slices, which results in extra terms like AH 2 
in equation Second, here we are dealing not with an 
initial but a boundary value problem, with associated 
boundary conditions for perturbations at the branes on 
the edges. After we derive the boundary conditions, we 
will calculate the KK spectrum of the eigenvalues m. 



III. BRANE EMBEDDING AND BOUNDARY 
CONDITIONS 

The embedding of each brane is described by w = 
w± + £±(x a ), where £± is the transverse displacement 
of the perturbed brane and w± is the position of the 
unperturbed brane. Holonomic basis vectors along the 

brane surface are = = (^ a , S^j , while the unit 

normal to the brane is ua = a^l + — £,,a5'Xj ■ The in- 
duced four-metric on the brane da 2 = q a bdx a dx b docs 
not feel the brane displacement (to linear order) and is 
conformally flat 



da 2 = a 2 {l- $) [-dt 



2 , e 2Ht 



dx 2 



(10) 



The junction conditions for the metric and the scalar field 
at the brane are 

dU 

[Kab - K,q a b] = U((p)q ab , [n-Vip] = — , (ff) 

a<p 

where the extrinsic curvature is defined by K a b — 
e fa) e (b) nA - B - We will only need its trace, which up to 
linear order in perturbations is 



2^-^. (12) 



For the background geometry (under the assumption 
of reflection symmetry across the branes), equations 
reduce to 



U iff 



U' 



a 2 T 6 ' a ± 2 ' 



(13) 



For the perturbed geometry, the traceless part of the 
extrinsic curvature must vanish in the absence of mat- 
ter perturbations on the brane. Since it contains sec- 
ond cross-derivatives of £, the brane displacement £ is 



severely restricted. Basically, this means that the oscilla- 
tory modes of brane displacement are not excited without 
matter support at the brane. While there could possibly 
be global displacements of the brane, they do not inter- 
est us, so in the following we set £ = 0. Of course, for 
the more complete problem which includes fluctuations 
&X °f the "inflaton" field on the brane, the displacement 
£ does not vanish. 

Using expression l|12|) for the trace of the extrinsic cur- 
vature, the first of equations (|ll|l gives us the junction 
condition for linearized perturbations at the two branes 
(a 2 <f>)' I = ±i U'a 3 5ip\ . However, this junction con- 

dition does not really place any further restrictions on 
the bulk field perturbations, as it identically follows from 
the bulk perturbation equations © and the background 
junction condition (|13|) . Rather, this junction condition 
would relate the brane displacement £ to the matter per- 
turbations on the brane if they were not absent. 

The second of equations l|ll|) gives us a physically rele- 
vant boundary condition for the bulk field perturbations 



(6<p'-<p'<f>)\ =±-U"aSw\ . 



(14) 



Using the bulk equations JSJ), this can be rewritten in a 
more suggestive form 



7 S(p 



3 m 2 
2 



AH 2 



a 2 $ 



atp ■ 



£VL — 42I zp a ij" 



(15) 

The eigenvalues m 2 of bulk perturbation equations sub- 
ject to the boundary condition i|15|) form a KK spectrum, 
which we find numerically. We considered several exam- 
ples of the potentials V and U± , and found no universal 
positive mass gap. Moreover, for the most interesting 
models we found negative m 2 . 

To understand the KK spectrum of m 2 , we make a 
simplification of the boundary condition l|15l) which will 
allow us to treat the eigenvalue problem analytically, and 
which well corresponds to a spirit of brane stabilization 
0- Indeed, rigid stabilization of branes is thought to be 
achieved by taking U" (i.e. the brane mass of the field) 
very large, so that the scalar field gets pinned down at the 
positions of the branes. In this case, the right hand side 
of (|15(l becomes very small, which leads to the boundary 
condition 



Sip\ 



0. 



(16) 



This by itself does not guarantee stability, or vanishing of 
the metric perturbations on the brane for that matter, 
as perturbations live in the bulk and only need to satisfy 
(|16fl on the branes. This poses an eigenvalue problem for 
the mass spectrum of the perturbation modes, which we 
study next. 



4 



IV. KK MASS SPECTRUM 

Unlike the situation with gravitational waves |22|. for 
the scalar perturbations there is no zero mode with 
to = 0, nor is there a "supersymmetric" factorized form 
of the "Schr6dinger"-like equation ©. To find the low- 
est mass eigenvalue, we have to use other ideas. Powerful 
methods for analyzing eigenvalue problems exist for nor- 
mal self-adjoint systems 25]. To use them, we transform 
our eigenvalue problem JSJ and 116|) into the self-adjoint 
form. While the second order differential equation 
is self-adjoint, the boundary conditions for u are not. 
Therefore, we introduce a new variable Y — u/z = a 2 $ 
and impose the boundary conditions (|16|l to obtain the 
boundary value problem 

VY = -(gYj + fY = XgY, (17a) 
Y'(w±)=0, (17b) 

where we have introduced the short-hand notation / = 
1/a, g = z 2 = (|a<y9' 2 ) , and A = m 2 + AH 2 . Since the 
boundary value problem i|17fl is self-adjoint, it is guar- 
anteed that the eigenvalues A are real and non-negative, 
A > 0. To estimate the lowest eigenvalue Ai of the eigen- 
value problem (TJJ, we apply the Rayleigh's formula |25j . 
which places a rigorous upper bound on Ai 



Ai < 



/ FVF dw 
J gF 2 dw 



(18) 



where F can be any function satisfying the boundary 
conditions (|17b(l . and does not have to be a solution of 
(|17afl . Taking a trial function F = 1, we have 



Jgdw' 



(19) 



This bound on the lowest mass eigenvalue is our main 
result: 



m 2 < -AH 2 



2 1^ 

3 f-H, 

J acp z 



(20) 



In practice, F = 1 is a pretty good guess for the lowest 
eigenfunction, so the bound iffijjl is usually close to satu- 
ration (up to a few percent accuracy in some cases) , as we 
have observed in direct computations using a numerical 
eigenvalue finder. 

The right hand side of equation (|20l) has the structure 
—AH 2 + m^(H), where the second term is a functional 
of H (including the implicit independence of the warp 
factor a) . In the limit of flat branes H — > we have only 
the second, positive term. In this limit our expression 
agrees with the estimation of the radion mass to 2 , fo r flat 
branes, obtained in various approximations [ill Il2l |20| . 
A non- vanishing H alters to 2 through both terms. The 
most drastic alteration of to 2 due to H comes from the 
big negative term —AH 2 . For the particular case of two 
de Sitter branes embedded in 5d AdS without a bulk 
scalar this negative term was noticed in |19|. 



V. TACHYONIC INSTABILITY OF THE 
RADION FOR INFLATING BRANES 

The most striking feature of the mass bound l|2(J|) is 
that to 2 for de Sitter branes is typically negative. Try- 
ing, for instance, to do Goldberger-Wise stabilization of 
braneworlds with inflating branes while taking bulk gra- 
dients ip 12 small enough to ignore their backreaction (as 
it is commonly done for flat branes) is a sure way to get a 
tachyonic radion mass: an estimate of the integrals gives 
to 2 < —AH 2 + 0((p' 2 ), which will go negative if the bulk 
scalar field is negligible (p' 2 <C H 2 . 

In what follows we consider two situations. In this sec- 
tion, we consider braneworld models where m 2 is negative 
and mostly defined by the first term —AH 2 in equation 
(|20|l . In the next section, we consider the case where 
both terms in equation H2(J|) are tuned to be comparable 
and the net radion mass is smaller than the Hubble pa- 
rameter | to 2 | < H 2 . In the last section we will discuss 
how these two cases may be dynamically connected. 

Suppose we start with a braneworld with curved de Sit- 
ter branes, and we find the mass squared of the ra- 
dion to be negative. The extra-dimensional eigenfunc- 
tion & m (w) is regular in the interval w_ < w < w+. Let 
us turn, however, to the four-dimensional eigenfunction 
Q m (t,x). Bearing in mind the evolution of the quan- 
tum fluctuations of the bulk field, we choose the posi- 
tive frequency vacuum-like initial conditions in the far 



past t -> -oo, f k (t) ~ ^2k elkn > V = I dte- Ht . For 
the tachyonic mode m 2 < the solution to equation Q 
with this initial condition is given in terms of Hankcl 
functions f^ m) (n) = ^-H \t)\ 3 / 2 H^ (kn), with the index 



fi = Y j + ^jpr ■ The late-time asymptotic of this solu- 
tion diverges exponentially as t — * oo (77 — * 0) 



fi m) (*) <* exp 




Ht 



(21) 



Thus the negative tachyon mass of the radion |to 2 | ~ AH 2 
leads to a strong exponential instability of scalar fluctu- 
ations $ oc e Ht . This instability is observed using a com- 
pletely different method in the accompanying BraneCode 
paper pj, where we give a fully non-linear numerical 
treatment of inflating branes which were initially set to be 
stationary by the potentials U± (</?) , and without any sim- 
plifications like approximating boundary condition l|15l) 
with JTHJ). 

Tachyonic instability of the radion for inflating branes 
means that, in general, braneworlds with inflation are 
hard to stabilize. From the point of view of 4d effective 
theory one would expect brane stabilization at energies 
lower than the mass of the flat brane radion TOg, which 
is roughly equal to the second term in (|2(J|I . If the energy 
scale of inflation H is larger than Too, H 2 ^> to 2 ,, this 
expectation is incorrect. 

Successful inflation (lasting more than requires 
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the radion mass m 2 to be not too negative 



as 



rn 



> 



IP 
20 ' 



da 2 = -(1 + 2$)rfP + (1 - 2^)e 2Ht dx 2 



(22) 



(23) 



This is possible if both terms in l|20|) are of the same 
order. In the popular braneworld models the radion mass 
in the low energy limit, mo, is of order of a TeV. For these 
models the scale of "stable" inflation would be the same 
order of magnitude, H ~ TeV. Although there is no 
evidence that this scale of inflation is too low, it is not a 
comfortable scale from the point of view of the theory of 
primordial perturbations from inflation. 

It is interesting to note that the system of curved 
branes may dynamically re-configure itself to reach a 
state where the condition (|22|) is satisfied. In the case 
of the bulk scalar field ip acting alone, for quadratic po- 
tentials U± suitable for brane stabilization, there may 
be two stationary warped geometry solutions Q with 
two different values of H . The solution with the larger 
Hubble parameter H might be dynamically unstable due 
to the tachyonic instability of the radion, which we de- 
scribed above. The second solution with the lower H 
which satisfies 11221) might be stable. A fully non-linear 
study of this model was performed numerically with the 
BraneCode and is reported in the accompanying paper 
0. It shows that, indeed, the tachyonic instability vi- 
olently re-configures the starting brane state with the 
larger H into the stable brane state with the lower H. 
This re-configuration of the brane system has a spirit of 
the Higgs mechanism. 

If we add an "inflaton" scalar field x located at the 
brane, its slow roll contributes to the decrease of H . 

Thus, for the "stable" brane we have a radion mass 
(|22|l . This condition includes the case when the radion is 
lighter than H, |m 2 | < H 2 . Even if the radion tachyonic 
instability is avoided, the light radion leads us to the 
other side of the story, a new mechanism of generation 
of scalar fluctuations from inflation associated with the 
radion. 



VI. INDUCED SCALAR METRIC 
PERTURBATIONS AT THE OBSERVABLE 
BRANE 

Suppose that the radion mass is smaller than H, 
\m 2 \ <C H 2 , so that from (|21|l we get the amplitude 

of the temporal mode function /£ (t) in the late time 
asymptotic frozen at the level (t) — J£ 3 / 2 ■ This is 
nothing but the familiar generation of inhomogeneities of 
a light scalar field from its quantum fluctuations during 
inflation. Therefore an observer at the observable brane 
will encounter long wavelength scalar metric fluctuations 
generated from braneworld inflation. 

The four dimensional metric describing scalar fluctua- 
tions around an inflating background is usually written 



where $ and \1/ are scalar metric fluctuations. The in- 
duced four-metric on the brane 111 Oil in our problem can 
be rewritten in this standard form 12: ill if we absorb the 
(constant) warp factor a(w+) in the redefined time t = at 
and spatial coordinates x = ax and rescale the Hubble 
parameter H = H/a. Then we see that the induced 
scalar perturbations on the brane are 



~ f 
2 



(24) 



The sign of the first equality here is opposite to what 
we usually have for 3 + 1 dimensional inflation with a 
scalar field. It implies that the 4d Weyl tensor of the 
induced metric vanishes, as the induced fluctuations are 
conformally flat. The conformal structure of fluctuations 
(|24|1 is typical j2^ for a R 2 inflation in the Starobinsky 
model |26|. It is not a surprise, because for the scale of 
inflation comparable to the mass mo of the flat brane 
radion we expect higher derivative corrections to the 4d 
effective gravity on the brane. Indeed, the massive radion 
corresponds to a higher derivative 4d gravity 20] . 

The amplitude and spectrum of induced fluctuations is 
defined by <&. From the mode decomposition 10 we get 



1VI4 



(25) 



where <& m (w+) is the amplitude of the extra-dimensional 
eigenmode at the observable brane, normalized in such a 
way that the fluctuations <I>(u>, i, x) are canonically quan- 
tized on the 4d slice, namely Mf J §^77 |$ m (w)| 2 dw = 1. 
The normalization M4 of the 4d mode functions follows 
from canonical quantization of the perturbed action J5J; 
the usual 4d Planck mass M p is expected to be recovered 
in the effective field theory on the observable brane [Tlj . 

The scalar metric fluctuations induced by the bulk 
scalar field fluctuations are scale-free and have the am- 
plitude fc 3 / 2 $fc cx -M-, with the numerical coefficient de- 
pending on the details of the warped geometry. The na- 
ture of these fluctuations is very different from those in 
(3 + l)-dimensional inflation, where the inflaton scalar 
field is time dependent. Induced scalar fluctuations do 
not require "slow-roll" properties of the potentials V 
and U±. The underlying background bulk scalar field 
has no time-dependence, but only y dependence. Thus, 
generation of induced scalar metric fluctuations from 
braneworld inflation is a new mechanism for producing 
cosmological inhomogeneities. 

If we add another, inflaton field \ localized at the 
brane, we should expect that fluctuations of both fields, 
the bulk scalar Sip and the inflaton 5x, contribute to the 
metric perturbations. We can conjecture that the net 
fluctuations will be similar to those derived in the com- 
bined model with R 2 gravity and a scalar field |27| . 
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VII. DISCUSSION 

Let us discuss the physical interpretation and the 
meaning of our result. Stabilization of flat branes is based 
on the balance between the gradient fa of the bulk scalar 
field and the brane potentials U(fa) which keeps <f> pinned 
down to its values fa at the branes. The interplay be- 
tween different forces becomes more delicate if the branes 
are curved. The warped configuration of curved branes 
has the lowest eigenvalue for scalar fluctuations around 
it 

m 2 = -AH 2 + ml(H) . (26) 

The term m§(i?) is a functional of H , and depends on the 
parameters of the model. If parameters are such that m 2 
becomes negative due to excessive curvature ~ H 2 , the 
brane configuration becomes unstable. This is analogous 
to an instability of a simple elastic mechanical system 
supported by the balance of opposite forces, which arises 
for a certain range of the underlying parameters. 

Tachyonic instability of curved branes has serious im- 
plications for the theory of inflation in braneworlds. It 
may be not so easy to have a realization of inflation in 
the braneworld picture without taking care of parame- 
ters of the model. Inflation where m 2 in l|26|) is negative 
and |m 2 | is larger than H 2 is a short-lived stage because 
of this instability. After inflation, the late time evolution 
should bring the brane configuration to (almost) flat sta- 
bilized branes in the low energy limit. This by itself re- 
quires fine tuning of the potentials V and U± to provide 
stabilization. Stabilization at the inflation energy scale 
requires extra fine tuning to get rid of the tachyonic ef- 
fect. 

Working with a single bulk scalar field, it is proba- 
bly not easy to simultaneously achieve stabilization not 
only at low energy, but also at the high energy scale of 
inflation, to insure that |m 2 | <C H 2 , and to provide a 
graceful exit from inflation. One may expect that intro- 
duction of another scalar field x on the brane can help 
to have stabilization both at the scale of inflation and 
in the low energy limit. If we can achieve brane sta- 
bilization during inflation by suppression of the tachy- 
onic instability, we encounter a byproduct effect. Light 
modes of radion fluctuations inevitably contribute to the 
induced scalar metric perturbations. Therefore the the- 
ory of braneworld inflation has an additional mechanism 



of generation of primordial cosmological perturbations. 
This new mechanism is different from that of the usual 
4d slow roll inflation. 

It appears that one of the most interesting potential 
applications of our effect is a mechanism for reducing the 
4d effective cosmological constant at the brane. Indeed, 
in terms of brane geometry, the 4d cosmological constant 
is related to the 4d curvature of the brane. Suppose 
we have two solutions of the background equations (J3J 
with higher and lower values of the curvature of de Sit- 
ter brane, which is proportional to H 2 . (The existence 
of two solutions for certain choices of parameters of the 
Goldberger-Wise type potentials used for brane stabiliza- 
tion can be demonstrated, see [llll0|.'l Suppose that the 
solution with the larger value of brane curvature is unsta- 
ble. Then the brane configuration will violently restruc- 
ture into the other static configuration, which is charac- 
terized by the lower value of brane curvature where the 
tachyonic instability is absent. The branes are flattening, 
which for a 4d observer means the lowering of the cos- 
mological constant. It will be interesting to investigate 
how this mechanism works for brane configurations with 
several scalar fields or potentials which can admit more 
than two static solutions. 

The problem of the cosmological constant from a 
braneworld perspective (as a flat brane) was discussed 
in the literature. There was a suggestion that the flat 
brane is a special solution of the bulk gravity/dilaton 
system with a single brane (2^, [3(j , the claim which was 
later dismissed 31]. In our setup, we consider two branes 
in order to screen the naked bulk singular ity, which was 
one of the factors spoiling the models [2(3, y(| • The new 
element which emerges from our study is the instability 
of the curved branes. 
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